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FOREWARD 


Interim  report  *  1  describes  progress  to  date  of  the  ELINT  Parameter¬ 
ization  and  Definition  Study.  The  introduction  describes  the  purpose  of  the 
study  and  explains  why  it  was  decided  to  model  a  superheterodyne  receiver  and 
defines  the  parameters  of  interest.  A  description  of  the  model  follows  in 
section  2.  A  flow  diagram  and  the  source  code  of  the  model  are  included  in 
sections  3  and  A  respectively.  The  report  is  concluded  with  documention  of 
results  and  conclusions  in  section  5. 


EXECUTIVE  SUMMARY 


This  Technical  Memorandum  was  prepared 
originally  as  part  of  the  Generic  ELINT/COMINT  Sensor 
Report  (FY-85/FY-86)  which  was  elmininated  under  the 
FY-87  statement  of  work  (SOW  #2) ,  undated  (delivered 
to  JPL  19  November  1986) . 

The  purpose  of  the  Generic  ELINT/COMINT  Sensor 
Report,  of  which  this  paper  was  intended  (in  its 
final  form)  to  become  part  of,  was  to  establish  a 
basic  superhetrodyne  receiver  based  sensor  model  and 
perform  simulations  with  it  to  determine  the  shaping 
or  coloring  of  the  statistical  distributions  of  the 
radar  free-space  signal  parametrics  by  a  typical 
sensor  prior  to  reaching  the  self-correlation 
processes.  It  was  also  intended  for  incorporation 
into  the  algorithm  test  bed  so  algorithms  could  be 
tested  with  realistic  distorted  data  rather  than 
unrealistic  stastically  pure  data. 

This  work  was  originated  in  support  of 
unanswered  questions  from  previous  self-correlation 
studies.  The  modeling  and  simulation  approach  was 
used  because  "live  date"  could  not  be  obtained. 

This  paper  is  being  published  because  it  was 
completed  in  FY-86  with  FY-86  funds  and  still  serves 
a  useful  function. 

There  were  no  unexpected  results  at  this  stage 
of  the  models  development  -  the  model  has  been 
verified  as  valid.  The  second  stage  of  this  study 
is  to  analyze  the  shaping  or  coloring  of  the  radar 
free-space  signal  parametrics  and  may  not  be 
completed  due  to  the  restrictions  of  the  current  JPL 
Task  Plan. 


/ 


ELINT  SENSOR  PARAMETERIZATION  AND  DEFINITION 
INTERIM  REPORT  *  1 


0.0  EXECUTIVE  SUMMARY 

1.0  INTRODUCTION 


2.0  RECEIVER  MODEL  DEVELOPMENT 

2. 1  SQUARE  LAW  DETECTOR  AND  LOW  PASS  DETECTOR 

2.2  VIDEO  AMPLIFIER 

2.3  CONSTANT  FALSE  ALARM  RATE  (CFAR) 

2.-4  PULSE  DETECTOR 

2.5  LIMITER 

2.6  DISCRIMINATOR 


3.0  RECEIVER  PROGRAM  FUNCTIONAL  FLOWCHART 

4.0  SOURCE  CODE 


5.0 


RESULTS  AND  CONCLUSIONS 


ELINT  SENSOR  PARAMETERIZATION  AND  DEFINITION 
INTERIM  REPORT  ♦  1 


1.0  INTRODUCTION 

The  purpose  of  this  report  is  to  define  the  statistics  of  the  observed 
parameters  of  interest  of  a  generic  ELINT  sensor.  The  approach  chosen  was  to  "build" 
a  computer  model  of  a  typical  ELINT  receiver  and  excite  it  with  input  signals  of 
known  distributions  and  then  measure  the  resulting  output  distributions. 

1.1  THE -SUPERHET* 

The  superheterodyne  receiver  has  many  advantages  to  offer  for  the  gathering  of 
electronic  intelligence.  It  is  the  most  widely  used  receiver  design  for  nearly  all  uses 
including  ELINT  applications.  A  basic  functional  block  diagram  is  shown  in  Figure  1.1. 
The  idea  is  to  use  a  local  oscillator  to  convert  the  incoming  signal  to  a  fixed 
intermediate  frequency  (IF)  by  the  mixing  process  called  heterodyning.  Thus  the  IF 
amplifier  need  operate  at  only  one  (lower)  frequency  and  it’s  operating 
characteristics  may  be  carefully  and  more  economically  controlled.  Bandwidths  for 
typical  narrowband  ELINT  receivers  are  approximately  20  MHz  for  Bpp  or  B|p  and  10 

MHz  for  By.  This  ensures  good  pulse  fidelity  for  nominal  radar  pulsewidths  of  I 

microsecond  or  greater.  The  simple  receiver  of  Figure  l.l  has  been  used  with  great 
sucess  since  World  War  II.  Recent  innovations  have  greatly  increased  the  basic 
capabilty  of  the  superhet.  By  rapidly  sweeping  the  local  oscillator  in  frequency  an 
amplitude  versus  frequency  display  may  be  generated  thus  creating  an  RF  spectrum 
analyzer.  Very  rapid  sweeping  creates  a  rapid  sweep  superhet  that  is  useful  for 
detecting  high  duty  cycle,  low  power  signals  in  the  presence  of  low  duty  cycle 
interfering  signals**.  In  addition  heterodyning  techniques  are  used  in  the 
implementation  of  compressive  and  microscan  receivers  and  many  instantaneous 
frequency  measurement  (IFM)  receiver  designs  utilize  superhet  frontends. 

So  we  see  that  the  study  of  superhet  receivers  is  a  good  place  to 
begin. 


"  R.  G.  WILEY,  ELECTRONIC  INTELLIGENCE:  THE  INTERCEPTION  OF  RADAR  SIGNALS 
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2.0  SENSOR  nOOEL  DEVELOPMENT 

Figure  2.1  is  a  functional  block  diagram  of  the  simple  receiver  to  be 

modeled.  We  first  assume  that  the  received  signal  is  a  narrowband  process  that 
is  statistically  unaltered  by  the  downconversion  to  the  IF  frequency,  Thus, 

the  signal  at  the  input  of  the  detector  and  limiter  is 

(I)  r(t)  =  s(t)  cosl^-u^t  ♦  <p(t)l  ♦  n(t)  where 

s(t)  =  transmitted  pulse  envelope  as  modified  by  the 
transmission  channel  characteristics 

<i>c  =  carrier  frequency 

<i)0  =  LO  frequency  assumed  to  be  noiseless 

o>if  =  o)c  -  <i)Q  =  IF  frequency 

<P(t)  =  signal  phase  at  output  of  IF  amplifier 
n(t)  =  additive  receiver  noise  plus  input  noise. 
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2.1  SQUARE  LAW  DETECTOR  AND  LOW  PASS  FILTER 


FIGURE  2.2  .  Functional  Block  Diagram  of  a  Square  Law  Detector 


Figure  2.2  depicts  functionally  a  square  law  detector.  Squaring  the  input 
( 1 )  above,  we  obtain 

r2(t)  =  s2(t)  cos2l2rtf  |Ft  ♦  <Kt)]  ♦  2s(t)  {cos(27Tf  |Ft  ♦  9(01)  n(t) 

*  n?(t) 

where  fJF  =  2Tt(<i)c-(i>0)  is  the  receiver  IF  frequency. 

Now,  using 

cos2x  =  1 1 2  (  1  ♦  cos  2x) 

(1.0  becomes 

r2(t)  =  1 1 2  s2(t)  ♦  1 1 2  s2(t)  cos  21  27tf  jFt  ♦  <p(t)l  ♦ 

2s(t)  {cost(2nf  ,F)t  ♦  <p(t)])  n(t)  ♦  n2(t) 
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Hence,  by  multiplying  the  received  signal  by  itself  we  have  obtained  '/2 

the  square  of  the  modulation  waveform  s(t)  plus  higher  frequency  terms 
consisting  of  cross  products  of  the  noise  and  the  signal  at  the  fund¬ 
amental  and  second  harmonic  frequencies  and  the  noise  term  nfy). 

Now,  by  low  pass  filtering  the  output  of  the  square  law  device  as  depicted  in 
the  above  figure  with  a  filter  sufficiently  wide  to  pass  s(t)  while  rejecting 
those  frequencies  above  the  lowest  component  of  s(t)  that  is  of  interest. 

One  obtains  at  the  output  of  the  detector  low  pass  filter  the  signal 

(2)  r v2(t)  =  V2  s*(t)  ♦  f?(t) 

where  n^t)  represents  the  low  frequency  components  of  the  noise  process 
n2(t). 


For  the  purpose  of  this  simulation  the  output  filter  of  the  detector  will 
determine  the  bandpass  response  of  the  video  channel.  Good  pulse  fidelity 
requires  that  the  filter  have  a  bandpass  adequate  to  pass  the  I2u>  harmonic 
with  less  than  75*  attenuation**.  Thus  the  filter  is  required  to  have  a  3db 
bandwidth  greater  than  or  equal  to  12  Mhz,  assuming  that  the  received 

pulse  width  is  greater  than  or  equal  to  Ipsec. 


**  R.  G.  WILEY,  ELECTRONIC  INTELLIGENCE:  THE  ANALYSIS  OF  RADAR  SIGNALS 
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2.2  VIDEO  AMPLIFIER 


The  video  amplifier  provides  additional  gain  and  bandpass  shaping  if  required, 
for  this  simulation  the  gain  will  be  set  to  K  with  an  infinite  frequency 
response.  Thus,  the  bandpass  response  of  this  portion  of  the  receiver  is 
determined  by  the  low  pass  filter  of  the  detector  and 

A(oj)  =  K 

where  K  is  selectable  as  an  user  input  into  the  simulation  and  defaults  to  the 
value  required  to  provide  unity  gain  to  the  detector  and,  thereby,  simplifying 
later  calculations. 


2.3  CONSTANT  FALSE  ALARM  RATE  (CFAR) 

The  CFAR  is  a  circuit  that  maintains  a  reference  voltage  level  proportional  to 
the  received  in  band  noise.  Thus,  it  may  be  used  to  vary  the  threshold  of  the 
pulse  detection  circuitry  in  such  a  way  to  prevent  false  alarms  due  to  varying 
noise  levels  at  the  input  of  the  receiver.  The  CFAR  constant  ,C,  is  adjusted  to 
provide  a  desired  false  alarm  rate(FAR).  The  results  reported  herein  were 
based  upon  a  calculated  FAR  =  10"4. 


2.4  PULSE  DETECTOR 

The  pulse  detector  circuitry  compares  the  output  of  the  video  amplifier  to  the 
threshold  determined  by  the  CFAR  circuit  and  strobes  a  pulse  present  signal  to 
the  parameter  estimation  circuit  when  the  video  signal  is  greater  than  the 
threshold.  A  threshold  is  chosen  to  provide  a  desired  probability  of  detection 
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2.5  LiniTER 


FIGURE  2.3.  FUNCTIONAL  BLOCK  DIAGRAM  OF  A  HARD  UNITER 


A  functional  block  diagram  of  a  hard  limiter  is  depicted  in  Figure  2.3.  The 
limiter  conditions  the  input  signal  to  the  frequency  discriminator  by 
amplifying  and  clipping  an  I.F.  output  signal  to  provide  essentially 
a  square  wave  input  to  the  discriminator  at  the  frequency  of  the  IF  signal.  If 
the  input  to  a  hard  limiter  is  x(t)  then  the  output,  y(t),  is 
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2.0  DISCRIMINATOR 


The  discriminator  provides  the  frequency  estimation  function  for  the  receiver. 
It  consists  of  a  zero  crossing  counter  that  counts  and  stores  the  number  of 
zero  crossings  during  the  pulse.  When  the  threshold  circuit  detects  a  pulse  it 
provides  a  strobe  that  initializes  and  starts  the  zero  crossing  counter.  At  the 
end  of  the  pulse  the  contents  of  the  counter  are  read  and  the  average  frequency 
of  the  l.F.  pulse  is  calculated  based  upon  the  estimated  pulsewidth. 
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r(k)  -  100r(k) 
If  r(k)  >-  0 
than  r(k)  - 1 
•1m  r(k)  -  -1 
tf  r(k)  >  r1(k) 
thifl  IFR  -  IFR  ♦  1 
rl(k)  -  r(k) 
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4.0  SOURCE  COOE  LISTING 


*6:34:07  PM 


♦THIS  PROGRAM  SIMULATES  A  SIMPLE  SUPERHETERODYNE  RECEIVER  THAT  UTILIZES  A 
SQUARE  LAW  DETECTION  PROCESS.  THE  IF  OUTPUT  SIGNAL  IS  FIRST  GENERATED  AS¬ 
SUMING  THAT  THE  RECEIVED  SIGNAL  IS  A  PULSED  RADAR  SIGNAL  OF  CONSTANT 
FREQUENCY,  RANDOM  PHASE,  AND  CORRUPTED  BY  AN  ADDITIVE  GUASSIAN  NOISE  CHAN¬ 
NEL. 


GENERATE  THE  RECEIVER  IF  SIGNAL 


K  -  SAMPLE  COUNTER 

KTOTAL  -  TOTAL  NUMBER  OF  SAMPLES  TO  BE  GENERATED 

SNR  -  SPECIFIED  SIGNAL-TO-NOISE  RATIO  OF  RECEIVED  SIGNAL 

R(K)  -  RECEIVED  SIGNAL  AT  OUTPUT  OF  IP.  AMPLIFIER 

FI  -  IP.  FREQUENCY 

DELT  -  THE  SAMPLE  INTERVAL  DELTA  T 

PHI  -  RANDOM  PHASE  OF  THE  RECEIVED  SIGNAL  R(K) 

S(K)  -  PULSE  ENVELOPE  OF  RECEIVED  SIGNAL 
PWID  -  PULSE  WIDTH 

M  -  NUMBER  OF  CFAR  INTEGRATION  SAMPLES 

C  -  CFAR  CONSTANT 

G  -  VIDEO  AMPLIFIER  GAIN 

IPP  -  PULSE  PRESENT  INDICATOR 

IFR  -  DISCRIMINATOR  ZERO  CROSSING  COUNTER 

CPRI  _  PRI  SAMPLE  COUNT 

IPW  -  PULSE  WIDTH  SAMPLE  COUNT 

PRI(N)  -  MEASURED  PULSE  INTERVAL  OF  NTH  PULSE 

PW(N)  -  MEASURED  PULSE  WIDTH  OF  NTH  PULSE 

FREQ(N)  -  MEASURED  FREQUENCY  OF  NTH  PULSE 


DECLARATIONS 


INTEGER  K,  IPP,  J,  I,  IFR,  CPRI,  IPW,  L,0,0 1  ,Q 
Integer**  TICKCOUNT,TICK  l,TICK2,toolbx 
REAL  Y(4200),X(4200),R(4200),  PRI (1000),  PW(1000),  FREQ' 1000) 
REAL  H(32),HL(24),RDET(4200) 

PARAMETER  (TICKCOUNT-Z'97550000)  !  toolbox  definitions 

DATA  (H(I),  I-U6)/-  5753*121E-02,  99027 196E-02,1.757335*5E-02, 


1-.65141 192E-02,  .13960525E-01,  22951469E-02, 

1  19994067E-0 1,7 1 3695&0E-02, 

1  -  39657363E-01, .  11260 114E-01,  .66233&43E-01, 

1-.  1 049722 3E-01,  65136133E-01,  -  12024993E-00, 

1-.29676577E00,304 109 17E00/ 

DATA  (HL(I),  I- 1, 12)/0.33?409 17E-02,0. 14936299E-0 1, 
10.10569630E-01, 

1 0.254 15067B-02,-0. 1 5929992E-0 1,  -0.34065343E-0 1, 

1-0.30112 177E-0 1,-0. 14629 169E-0 1, 

1  0.4006954 1E-0 1.0. 1 15407 13E00.0. 13650752EOO, 

IO.233546O6HOO/ 

PRINT  V  PULSE  PULSEWIDTH  PRI  FREQUENCY 
1  SNR' 

PRINT  V  PREDETEC 

1TION  POSTDETECTION' 

OPEN(UNIT«16JILE  «  ’RCVDATA*) 

DO  111  1-17,32 

L-33-I 

H(I)-H(L) 

1  CONTINUE 
DO  101  1-13,24 
L-25-1 

HL(I)-HL(L) 

1  CONTINUE 

WR1TE(9  *)  ENTER  NUMBER  OF  PULSES,  MINIMUM  SNR,  MAXIMUM  SNR 

COUNTERS  AND  FLAGS****************** 

READ(9  ,*)  O.SNRMIN.SNR 
TICK1  -  toolbx(T I CKCOUNT ) 

DELT  »  2.5E-06 
FI  -  10E06 
KTOTAL  -  4000 
PWID  -20E-06 
p  01-0 
Q-0. 

SNSUM  -  0. 

SNSUM1  -0. 

PM-0. 

II  -  1 

I  IPP  «  0 
IPW  -  0 

II -1 

IFR-0 
PWRSIG  -  0. 


p  i.»  r  f 


>.v 

ft?* 
> 


mm 


’I  *  I,1  «,*  M 


y.y.v:  v.  y.v.v.v.v.  .v.v.v..-. 


jEgSl 

5* 


V» 
O.S> 


IwwfwSwwwl?* 


PWRNOS  •  0. 

PRESIG  -  0 
PRENOS  ■  0. 

KF-5555 

*  WRITE<9,*)  -ENTER  SAMPLE  INTERVALS  F.  FREQUENCY  DUMBER  OF 

*  1  SAMPLES1 

*  READ(9 i  *)  DELTp  l^TOTAL 
PI  -3  141593 

♦GENERATE  RANDOM  PHASE  PHI 
SIGMA  -  1  ./SQRT(SNR) 

PHI  -  RAND(O)  !  SEED  THE  RANDOM  NUMBER  GENERATOR 

PHI  *  (RAND(l)  -  5)*2  0*PI 
IDELAY  -  2000 
L  «  INT(PWID/DELT) 

DO  1 1  K«  l,KTOTAL  ♦  130 

TN  *  RAND(l)-  5 

XN  -  SQRT(LOG(  l./(TN*TN))) 

XN  -  XN  -  (2.30753  ♦  2706l*XN)/(l.  ♦  99229*XN  ♦  044dl*XN*XN) 

XN  «  XN*SIGMA 
1FCTN.LT.0)  XN«-XN 
•GENERATE  PULSED  IF  SIGNAL 

IF  (K  LE.  IDELAY  OR.  K  GE  L  ♦  IDELAY)  THEN 

R(K)  -  XN 

ELSE 

R(K)  ■  1.41 42*COS(2*PI*F  1*K*DELT  ♦  PHI)  ♦  XN 
END  IF 

1 1  CONTINUE 

*  BANDLIMITED  ZERO  MEAN  WHITE  GAUSSIAN  NOISE 

*  BANDPASS  NOISE  TO  ♦/-  F 1/2 

*  FINITE  IMPULSE  RESPONSE  (FIR)  LINEAR  PHASE  BANDPASS  DIGITAL  FILTER 

*  -  10DB  POINTS  ARE  14  AND  39  OF  THE  NYQUIST  FREQUENCY 

DO  102  K-32XTOTAL+130 
TOO  •  0 
DO  102  J  -  1,32 

102  YOO  -  YOO  ♦  H(J)*R(K-]+l)  ©ANDPASS  LIMITED  PULSE  IF  PLUS  NOISE 
D0  21K-1,KT0TAL*99 

Y(K)  -  Y(K>3 1  )*Y(K*3 1 )  (PERFORM  SQUARE  LAW  DETECTION 

IF(K.LE.IDELAY-32  OR.  K.GE.IDELAY*L)  THEN 
PRENOS  -  PRENOS  ♦  YOO*Y(K)/(KTOTAL  ♦  67  -  L) 

ELSE 

PRESIG  -  PRESIG  ♦  T(K)*Y(K)/L 
END  IF 

2 1  CONTINUE 


RDET(l) 


PRESNR  ■  PRESIG/PRENOS 
PRESNR  -  PRESNR 
PSNRDB  -  10*LOG10(PRESNR) 

****************LOW  PASS  FILTER  DETECTOR  OUTPUT******** 
***********-10  DB  POINT  IS  .05  OF  THE  NYQUIST  FREQUENCY** 
RDET(l)  -  0. 

DO  202  K*  1,100 

RDET(l)  »  RDET(l)  +  01*Y(K) 

CONTINUE 

DO  203  K-2,KT0TAL 

RDET(K)  -  RDET(K-l)  ♦  01*(Y(K*99)-Y(K-1)) 

IF(K.LE.IDELAY- 130  OR.  K.GE.IDELAY+L)  THEN 
PWRNOS  -  PWRNOS  ♦  RDET(K)*RDET(K)/(KTOTAL  -  L- 130) 
ELSE  IF  (K.GE.IDELAY-130  AND.  K.LE.IDELAY+L)  THEN 
PWRSIG  -  PWRSIG  ♦  RDET(K)*RDET(K)/(L- 130) 

END  IF 
CONTINUE 

SIGNOS  «  PWRSIG/PWRNOS 
SIGNOS  «  SIGNOS  -  1.0 
SIGNOSDB  *  10  *L0G  lO(SIGNOS) 

SNSUM  -  SNSUM  ♦  SIGNOS/O 
SNSUM 1  -  SNSUM  1  ♦  PRESNR/0 


T  *  SQRT(16.4*PWRNOS) 

DO  400  K-l^TOTAL 

IF(  RDET(K)  LT.  T  AND.  1PP-0  AND.  K-KTOTAL)  GO  TO  206 
IF(  RDET(K)  GE.  T  AND.  IPP-0  )  THEN 
KF-K 

T-T/2.  !  FIRST  SAMPLE  OF  NEW  PULSE 

IPP-1 

PRIG  1)  -  (2000.  ♦  K)*DELT*  1000000. 

IPW  «  IPW  ♦  1 

ELSE  IF(  RDET(K)  GE.  T  AND.  IPP- 1  )  THEN 
IPW  -  IPW  ♦  1 


ELSE  IF(  RDET(K)  LT.  T  AND  IPP- 1  )  THEN 
PW(1 1)  -  IPW*DELT*  1000000. 

GOTO  401 
END  IF 

400  CONTINUE 

401  1F(KFNE.5555)  THEN 

DO  204  K-KF,KF*IPW 
IF(R(K).LT  0)  THEN 
ROO--1. 


ELSE 
ROO  *  1. 

END  IP 

204  CONTINUE 
END  IF 

DO  205  K-KF^+IPW 

205  IF(R(K)NER(K+D)  IFR  *  IFR  +  1 
FREQ(Il)  -  IFR/PW(I l)/2. 

*  PRINT  5002  ,1 U>W(I  DPRKI  OjREQd  DJPSNRDB,  SIGNOSDB 
*5002  FORMAT(  1X,I 4,3X,F  1 0.4,3X,F  1 0.4,6X,F  1 0.4, 4X,F  1 0.4.2XJ  10.4) 

GO  TO  207 

206  PM  =  PM  ♦  l./O 

*  PRINT  5003, 1 1 XSIGNOSDB 

*5003FORMAT(1X,I4,6X,16HMISSED  DETECTION,  5X,  1 2HTHRESHOLD  =  J10.4, 

*  15X,6HSNR  =  ^10.4) 

207  II  -  II ♦  1 
01-01-1 

IF  (Ol.GT.O)  GO  TO  99 
SNSUM  =  10*LOG  lO(SNSUM) 

SNSUM1  *  10*LOG10(SNSUM1) 

WRITE(6  *)  AVERAGE  PREDETECTION  SNR  -  ,SNSUM  1/ 

1AVERAGE  POSTDETECTION  SNR  -  \SNSUM 
PD  *(1.  -  PM)*100. 

WRITE(6*)  PROBABILTY  OF  DETECTION  «  J>D, '  PER  CENT' 

TICK2  •  toolbxflTCKCOUNT ) 

TIME  «  FLOAT(TICK2  -  TICK  l)/60. 

WRITE(6  *)  ELAPSED  TIME  -  '/TIME/  SECONDS' 

SNR  *  SNR  -.10 

IF(SNR-SNRMIN  >  0.)  GO  TO  97 

*  D0  91K-UTT0TAL 

*  X(K)  -  RDET(K) 

*  Y(K)  -  0.0 

*91  CONTINUE 

*  CALL  FFT4(X,Y,4096,6) 

*  NTOTAL  -  KTOTAL/2 

*  DO  40  K»  1  JfTOTAL 

*  RDET(K)  -  X(K)*X(K)  ♦  Y(K)*Y(K) 

*  RDET(K)  «  10  *LOG  10(RDET(K)) 

*40  CONTINUE 

**************************  AyERAGE  ARRAY  ********************************* 

*  LTOTAL  -  NTOTAL/4 

*  DO  50  K-l.LTOTAL 

*  Y(K)  •  0.0 


20 


50  CONTINUE 
L«1 
M»4 

DO  60  K«  1,LT0TAL 
DO  70  J-L,M 

1  Y(K)  -  Y(K)  ♦  RDET(J)/4. 

70  CONTINUE 
L«L+4 
M-M+4 

60  CONTINUE 

[  WRITE(  16,50 10)  <RDET(K),K-  l^TOTAL) 
!5010FORMAT(2XJ154) 

WRITE(9 i  *)  ‘ENTER  1  TO  CONTINUE1 

READ(9,*)Q 

IF(Q«1)  GO  TO  96 

STOP 

END 


FUNCTION  Rand(IX) 

* 

*  This  random  number  generator  is  a  variation  on  Tausworthe 

*  generator  described  in  “Solution  of  Statistical  Distribution 

*  Problems'  by  H.  0.  Hartley  in  Statistical  Methods  for  Digital 

*  Computers  Voi.  Ill,  edited  by  Enstein,  Ralston,  and  Will 

*  (John  Wiley  and  Sons  1977) .  The  only  modification  of  consequence 

*  is  the  ability  to  reseed  the  generator  with  the  system  tickcount. 

* 

*  The  function  returns  a  real  value  between  0  and  1. 

*  Rand(O)  will  reseed  the  random  sequence  using  system  tickcount. 

*  Rand(l)  will  use  previous  calls  values  as  seeds  for  next  number 

*  in  the  sequence.  Always  using  Rand(  1 )  will  generate  a  specific 

*  random  sequence  based  upon  starting  values  internal  to 

*  the  function. 

* 

Real*4  Factor  Jfu  l,Ku2 

Integer* 4  toolbx,  B1TX0R,  TICKCOUNT  !  toolbox  definitions 
Integer*4  KU3.LU  1,KC,N2TM,N2TCM 
Equivalence  (KU2,KU3) ,  (KU  1,LU  1) 

PARAMETER  (BITXOR-Z‘65992000  )  !  toolbox  definitions 

PARAMETER  (TICKCOUNT-Z,97560000)  !  toolbox  definitions 

*  Note  that  the  parameter  statement  assigns  a  value  to 

*  the  indicated  constant  label.  In  this  case,  it  is  the  indey 


*  into  a  table  of  trap  addresses  contained  in  TOOLBX.SUB  which 

*  converts  the  call  to  a  trap  call.  A  list  of  these  assignments 

*  is  given  in  the  file  TOOLBX.PAR. 

* 

*  You  must  explicitly  tell  MacFortran  to  save  the  values  of 

*  local  variables  across  sucessive  calls  of  the  subroutines. 

* 

SAVE 

Data  Factor/0.46566l3e-9/ 

Data  N2TCM/Z000400007 

Data  KU  UC^2TM/z*40000003  ,z7fffffff ^  00002  0007 

*  KU  1  enters  with  U(I)  uniform  random  variable 

*  KU2  leaves  with  U(M)  uniform  random  variable 

*  KCU1  complement  of  KU  1 

*  KCU2  complement  of  KU  2 

*  KC  complementing  constant 

*  LXXXO  logical  equivalences  of  above  LXXX 

*  N2TM  2**M  where  Mis  shift  factor 

*  n2TCM  2**P-M  where  P  is  word  size 

*  Factor  float  value  of  2**P 

*  XOR  is  the  exclusive  or  operator 


if  (IX-0)  KUl«tOOlbx(TICKCOUNT) 

KU  l-toolbxfBITXORJUJ  1,LU  1/N2TM) 

KU2  -toolbx(BITXOR,KU  1,(LU  1*N2TCM)  and.  KC) 

Rand«FLOAT(KU3)*Factor 

KU1-KU2 

RETURN 

END 


C  A  COOLEY -TUKEY  RADIX-4  DIF  FFT  PROGRAM 
C  COMPLEX  INPUT  DATA  IN  ARRAYS  X  AND  Y 
C 
C 
C 

C . 


SUBROUTINE  FFT 4  (X,Y,N,M) 
REAL  X(  1024),  Y(  1024) 
REALM  E 


C 


N2  -N 


MAIN  FFT  LOOPS 


N1  -N2 
N2  -  N2/4 

E-6.2B31 65307 17953&/N1 
A  ■  0 

DO  20  J  -  1,N2 
B  ■  A  +  A 
C  *  A  +  B 
C01  -COS(A) 

C02  *  C0S(B) 

C03  «  COS(C) 

SI  1  *  SIN(A) 

SI  2  *  SIH(B) 

SI  3  -  SIN(C) 

A  « J*E 
DO  30  I*J,  N,  N1 
1 1  •  I  ♦  N2 

12  -I1+N2 

13  -  12  +  N2 

R1  »  X(I)  ♦  X(I2) 
R3«X(I)-X(I2) 

51  ■  Y(I)  +  Y(I2) 

S3  *  Y(I)  -  Y(I2) 

R2  »  X(1 1)  +  X(I3) 

R4  ■  X(1 1)  -  X(1 3) 

52  ■  Y(1 1)  ♦  Y(I3) 
S4-Y(I1)-Y(I3) 

X(I)  «  R1  ♦  R2 

R2  -R1  -R2 
R1  -R3-S4 
R3-R3  +  S4 
Y(I)  -  SI  ♦  S2 

52  *  SI  -  S2 
SI  »  S3  ♦  R4 

53  •  S3  -  R4 

X(1 1)  ■  C01*R3  ♦  SI  1*S3 
Y(1 1)  •  C0I*S3  -  SI  1*R3 
X(I2)  ■  C02*R2  ♦  SI2*S2 
Y(I2)  *  C02*S2  -  SI2*R2 
X(1 3)  -  C03*R  1  ♦  SI  3*S  1 
Y(I3)»C03*S1  -SI3*R1 
CONTINUE 
CONTINUE 


CONTINUE 

. DIGIT  REVERSE  COUNTER - 

J-l 

N1  -N  -  1 
DO  104  I- 1,  N 1 
IF  (I.GE.J)  GOTO  101 
R1  «X(J) 

X(J)*X(I) 

X(I)  a  R1 
R1  «Y(J) 

Y(J)  -  Y(I) 

Y(I)«R1 

K-N/4 

IF  <K*3  GE.J)  GOTO  103 
J.J-K*3 
K-K/4 
GOTO  102 
J-J  +  K 
CONTINUE 
RETURN 
END 
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S.O  RESULTS  AND  CONCLUSIONS 


The  previous  pages  of  this  report  have  documented  the  development  of  the 
Sensor  model.  In  this  section  model  validation  and  verification  will  be 
demonstrated.  In  order  to  validate  the  model  a  received  signal  disturbed  only  by 
Guassian  noise  was  assumed  and  the  signal  processing  was  examined  at 
intermediate  points  within  the  model  to  see  if  the  model  was  behaving  as  one 
would  expect. 

The  assumed  signal  parameters  for  the  validation  runs  were  Fjp  =  to  MHz  and  a 
pulsewidths  of  10,20,  and  50  jisec. 

A  FFT  routine  was  used  to  observe  the  output  of  the  IF  bandpass  filter.The 
results  of  two  runs  of  SNRs  of  0  db  and  10  db  (SNR  measured  at  input  of  IF  filter) 
are  depicted  in  Figure  5.1.  As  expected  the  output  is  a  bandpass  spectrum  centered 
at  10MHz  with  a  spectral  width  of  10  MHz.  Next  the  output  of  the  detector  lowpass 
filter  was  observed.  This  is  illustrated  in  Figure  5.2.  Again  we  see  what  one  would 
hope  to.  A  lowpass  spectrum  of  spectral  width  of  about  5  Mhz  (adequate  to 
accurately  reproduce  pulsewidths  of  2  jisec  and  greater)  and  the  signal  energy 
concentrated  at  the  lower  frequencies. 

Figures  5.3,  5.4,  and  5.5  depict  time  domain  plots  of  the  detector  output. 

The  received  pulse  SNR  was  20  db,  13  db,  and  3  db  respectively  in  the  IF  bandwidth. 
The  stronger  pulses  were  detected  and  the  pulsewidth  of  50  jisec  correctly 
measured.  At  a  SNR  =  3db  false  detection  occurred  and  pulsewidth  was  incorrectly 
measured.  The  plots  as  shown  have  been  smoothed  to  facilitate  the  plotting 
process.  This  enhances  the  SNR  as  seen  in  the  plots  and  the  detector  was  working 
at  a  SNR  of  about  I3db  less  than  one  would  estimate  by  looking  at  the  plots  alone. 


5.1.  3POTHM.  PUTTS  or  THE  SIMM.  AT  TNI  OUTPUT  OF  TNI  I.F.  FILTER 
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DETECTOR  FILTER  OUTPUT 


0\  0.  0.  0.  1.  I.  I.  2.  2.  2.  3.  3.  3.  4  4  4 

-lol  "^3  6  9  2  5  8  1581  47036 

1  \2  35678  1235678 

2  5  7  257  257  257 

““T  5  \  5T\  5  5  5  5  5  5 


RELATIVE 


POWER  DB 


12  3  5  6  7  8 
2  5  7  2  5  7 
5  5  5  5 


FREQUENCY  MHZ 
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FIGURE  5.3.  DETECTED  PULSE  WITH  MGH  SNR 
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FIGURE  5.4.  DETECTED  PULSE  OF  LOWER  SNR. 
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The  model  was  exercised  many  limes  with  received  signals  of  varying  SNRs  and 
then  the  measured  probabilty  of  detection  was  compared  to  theoretical  as  docu¬ 
mented  in  many  standard  references*.  A  plot  of  SNR  vs.  Probability  of  Detection  for 
this  model  is  presented  in  Figure  Eight. 

nUMM.TY  OF  KTECTRM  V8. 9Nt 
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FBURESI.  PR0MM.TY  OF  DCTECTION  V8.  SNR  PLOT  FOR  TK  EL»IT  SOBOR  MOOa 

As  can  be  seen  performance  is  comparable  at  high  SNRs  but  degrades  rapidly  as 
SNR  decreases.  This  is  because  the  detection  process  modeled  is  not  matched  filter 
detection  and  is  not  optimum  at  lower  SNRs  but  is  instead  designed  to  reproduce  a 
detected  pulse  with  sufficient  fidelity  to  provide  good  pulse  parameter  measure¬ 
ment  capabilty. 

In  conclusion,  then,  it  can  be  said  that  the  sensor  model  at  its  present  state  of 
development  performs  as  expected  when  excited  with  signals  plus  Gaussian  noise 
and  will  provide  a  tool  by  which  the  output  distributions  of  other  types  of  signals 
may  be  studied. 
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SEE  FOR  INSTANCE  “RADAR  TARGET  DETECTION  HAND  BOOK  OF  THEORY  AND 
PRACTICE  m  BY  MEYER  AND  MA  YER 
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